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Sensing under Quantization 
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Abstract 

In cognitive radio, the cooperative spectrum sensing (CSS) plays a key role in determining the 
performance of secondary networks. However, there have not been feasible approaches that can analyt- 
ically calculate the performance of CSS with regard to the multi-level quantization. In this paper, we 
not only show the cooperative false alarm probability and cooperative detection probability impacted 
by quantization, but also formulate them by two closed form expressions. These two expressions enable 
the calculation of cooperative false alarm probability and cooperative detection probability tractable 
efficiently, and provide a feasible approach for optimization of sensing performance. Additionally, to 
facilitate this calculation, we derive Normal approximation for evaluating the sensing performance con- 
veniently. Furthermore, two optimization methods are proposed to achieve the high sensing performance 
under quantization. 

Index Terms 

cooperative spectrum sensing, cognitive radio, fusion rule, quantization, optimization. 

I. Introduction 

Recently, extensive studies found that the scarce spectrum resource licensed to primary users 
(PUs) are under-utilized significantly [l]-[4]. To alleviate this issue, cognitive radio (CR) technol- 
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ogy has been brought as a promising method to make efficient use of the limited radio spectrum. 
In CR, the secondary users (SUs) are permitted to dynamically utilize the idle frequency spectrum 
by sensing the reasonable access opportunity. Thus, the accuracy of spectrum sensing plays a 
key role in determining the performance of CR system. To identify the access opportunity with 
high correctness, the cooperative spectrum sensing (CSS) has been widely concentrated since it 
greatly improves the sensing performance. 

A. Related Works 

For CSS, the fusion center makes a final sensing decision following a certain fusion rule, 
based on the total sensing information reported from each involved cooperative user. Commonly, 
the CSS with fusion rule are generalized into two major categories: i) pre-processing based 
fusion rule [5]— [13]. The fusion center processes the measurements (or test statistics), which are 
generated by cooperative users pre-processing the observation, to make final sensing judgment; 
ii) direct-forward based fusion rule [14]— [16]. The fusion center processes the total received 
samples which are forwarded from each cooperative user to make a final sensing decision. The 
direct-forward based fusion rule requires a large number of overhead for cooperative users to 
feed back their collected samples, which means the gain from cooperation can be exhausted by 
the overhead on data fusion. Hence, when compared with the direct-forward based fusion rule, 
the pre-processing based fusion rule is more valuable in CR. 

As shown in [5], the pre-processing based fusion rule is called as the multi-bit fusion rule based 
on multi-level quantization. In digital communication, it is not avoidable that the measurements 
(or test statistics) are quantized and converted into a binary output code for transmission. 
Hence, the test statistics is transmitted by multiple bits and in quantized manner all the time. 
On the other hand, both the accuracy and the temporal response are crucial in CSS. The 
temporal response depicts how fast the primary user can be detected out since the presence 
of primary user, and it is mainly determined by the quantity of overhead spent on reporting 
the test statistics. As pointed out in [17], the cooperation overhead makes CSS impractical 
under large number of cooperative users. Hence, the less overhead on reporting the test statistics 
achieves the better temporal response. This is the reason why there is an amount of researches 
concentrating on the CSS with counting rule (or called as decision fusion rule), which is 
also interpreted as the bi-level quantization based fusion rule [18]. Briefly, the number of 
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quantization levels is as small as possible from the perspective of reducing the cooperative 
overhead. However, the limited quantization existing ubiquitously in wireless communication 
introduces quantization error deteriorating the performance of CSS. Thus, it is necessary to 
study the impact of quantization on the performance of CSS. 

For the CSS with multi-level quantization, current researches have not provide a feasible 
approach to ensure the false alarm probability and detection probability of CSS can be tractable 
by closed form expressions [5], [14]. Although [5] showed closed form functions to evaluate the 
final sensing performance (including cooperative false alarm probability and cooperative detection 
probability), there are unknown variables in these two functions of which the values can only be 
obtained by exclusive search. Thus, the computation complexity on evaluating the performance of 
final decision in [5] is very difficult, causing the related performance evaluation improper online. 
[14] indicated a two-bit fusion rule based CSS scheme where each cooperative user uses three 
thresholds to divide the value range of the observed test statistics into four regions with different 
weights and employs two-bit report to express the region where the value of test statistics falling. 
However, in [14], the weights are preferred unreasonably and the decision threshold cannot be 
adjusted. As presented in [14], the fusion rule with multi-level quantization is left for the future 
work. Since the performance of CSS with quantization error cannot be computed efficiently and 
presented analytically, the current CSS schemes related to multi-level quantization cannot apply 
the Neyman-Pearson (NP) criterion and Constant Detection Rate (CDR) criterion reasonably, 
and the optimization of CSS performance cannot be performed efficiently. 

B. Contributions 

In CR, the common CSS's application scenario, called as homogeneity scenario in this paper, 
is that the distances between any SUs are small compared with the distance from any SU to the 
primary transmitter. Thus, the current literatures [6]— [13] regard the averaged received signal to 
noise ratio (SNR) at each SU are identical. This paper studies the performance of CSS with 
multi-level quantization in homogeneity scenario. 

The quantization theory is employed to have an insight into the current multi-bit fusion rule. 
Then the multi-bit fusion rule is interpreted as the multi-level quantization based fusion rule. 
From the perspective of quantization theory, the study of final sensing performance is to analyze 
the impact of quantizing the test statistics on the data fusion. Based on this opinion, we derive 
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the cooperative false alarm probability and cooperative detection probability given by closed 
form expressions. Especially, this paper does not specify a certain type of detector for the 
derivation of expressions. Consequently, the derived expressions have the following merits: /) 
they are eligible for arbitrary number of quantization levels; ii) they are universal for different 
types of detector; Hi) they guarantee both the cooperative false alarm probability and cooperative 
detection probability can be tractable and optimized reasonably; iv) moreover, the derived closed 
form expressions can be used to analyze the impact of quantization error in different quantization 
scheme, for example, the uniform quantization scheme. 

To facilitate the related performance evaluation, the discrete Fourier transform (DFT) method 
is introduced to simplify the derived two closed form expression and reduce the corresponding 
computation complexity. Additionally, the Normal approximations for the obtained two closed 
form expressions are derived to estimate the CSS sensing performance conveniently. Based on 
the derived Normal approximations, the optimization of sensing performance are programmed 
to a linear combination problem which can be solved by semidefinite programs efficiently [16]. 
Furthermore, when considering the Lloyd [19] quantization scheme, an optimization method is 
presented by using quasi-convex theory. 

C. Paper Structure 

The rest of this paper is organized as follows: Section II describes the system model, and 
formulates the fusion problem in current literatures. Sequentially, Section III presents the per- 
formance evaluation of CSS with multi-bit fusion rule. In Section IV, the quantization schemes 
and related optimization methods are discussed. Followed by conclusion, numerical studies are 
indicated in Section V. 

II. Preliminary and Sensing Model 

A. Sensing Model 

In CR, SUs are required to sense the presence of distant primary transmitter to ensure the QoS 
of primary link. Considering a secondary base station and related accessing SUs cooperatively 
sense the presence of primary user, the CSS problem is modeled by Fig. 1. The sensing objective 
is to properly determine one of two hypotheses H (denote primary user off) and Hi (denote 
primary user on) by analyzing the received signal X u where % e {1, 2, • • • , K} denotes the z'-th 
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SU. The K SUs pre-process the observed signal to generate the related test statistics of which 
the probability density functions (pdfs) are f(x\H ) and f(x\Hi) for H and Hi respectively. 1 
According to the value of test statistics T; online, each cooperative SU generates its own report 
Ri, and then they feed back R { to secondary base station. After the collection of reports Ri, the 
secondary base station as the fusion center utilizes the received Ri to make a final decision R Q 
for the whole cooperative secondary cluster. 

In this paper, consider the received signal Xi is independent and identically distributed (i.i.d.) 
for different cooperative SU, and f^°° f(x\H )dx < f(x\Hi)dx for Ve G (—00, +00). 
These requirements hold for common detectors employed in CR. 

B. Preliminary 

For conventional-decision-fusion-rule based CSS, the report Ri of i-th SU is generated as 
following (1), 

fo, if Ti < e 

Ri(Ti) = I (1) 
[l, if Ti > e 

where e is the sensing threshold, and Ti is the test statistics at the z-th SU. And the fusion 
center makes final sensing decision according to the number of SU reporting the occurrence of 
primary user. In fusion center, the final sensing decision strategy R (-) is given by 

(h , if Y:tiRi< N 
R (Ri) = < (2) 

(Hi, if E?=iRi>N 
where jV is termed as decision threshold and N e {1, 2, • • • , K}. 

According to the probability theory, the related cooperative false alarm and detection proba- 
bilities under decision fusion rule are given by 

P(T > e \H ) k (l - P(T > e \H )) K - k , (3) 
P(T > e |#i) fc (l - P{T > eol^i))^"" (4) 

'in statistical hypothesis testing, a hypothesis test is typically specified in terms of a test statistic, which is a function of the 
sample [20]. 
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where Pf and P d denote the cooperative false alarm probability and cooperative detection 
probability respectively; T is short for Ti\ P(T > eo\H ) and P(T > e \Hi) are the false 
alarm probability and detection probability of individual cooperative user respectively. 

C. LRT based Multi-bit Fusion Rule 

For CSS in homogeneity scenario, [5] presents a certain approach to calculate false alarm 

probability and detection probability regarding multi-level quantization. Let I denote the number 

of levels in quantization. The whole range of test statistics is separated into / regions: region 

0, region 1, • • • , region / — 1. In addition, {k , ki, • • • , ki-±} mean the number of SUs whose 

test statistics falls in region 0, region 1, ■ • ■ , region I — 1 respectively. Thus, according to [5], 

when the likelihood ratio test (LRT) is used at the fusion center and S represents the set of all 

{k ,ki,--- , ki-i} satisfying that Y^i=o ^ = ^ an< ^ (^), the final probability of detection and 

probability of false alarm are given by (6) and (7) respectively. 

P(T < eo|ffi) fe °P(e < T < e^H^ ■ ■ ■ P(q_ 2 < TjH^ C f P(H ) 
P(T < e \H ) k oP(e <T< e^H^ ■ ■ ■ P(e,_ 2 < T\H ) k ^ > C m P(fTi) 

where Cj and C m are weights. 

{koM,- M-i}es , (6) 

xP(e < T < e^H^ ■ ■ ■ P(e,_ 2 < T^f^ 

p >= E Wh s.:. kl _^ T <^ 

{k M,- M-i}eS . (7) 

xP(e < T < exl^o)" 1 • • • P(Q-2 < Pl^o)^ 1 
For (6)(7), the set S is obtained by substituting the every permutation of {ho, fci, • • • , h-i} 
into (5). Evidently, it cannot be efficient to acquire S for inequality (5). Hence, the computation 
complexity of final detection probability and false alarm probability is high, when regarding the 
computation of acquiring S. Additionally, the analysis of the sensing performance (Pf,P d ) and 
the choice of optimal threshold controlling (Pf,P d ) are mathematically intractable in current 
literatures generally. 
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III. Performance of CSS 



This paper treats the conventional multi-bit fusion rule as the multi-level quantization based 
fusion rule from the perspective of quantization theory, and then studies the performance of CSS 
with multi-level quantization mathematically. To facilitate the understanding of our results, the 
decision fusion rule based CSS is re-interpreted firstly as an instance of using quantization theory, 
when compared with the conventional interpretation in Section II-B. Sequentially, the cooperative 
false alarm probability and cooperative detection probability of the multi-level quantization are 
derived and explained in detail. 

A. Interpretation of Decision Fusion Rule from Quantization Theory 

In the view of quantization theory, the decision generated in SU can be regarded as a special 
result of multi-level quantization where 1 — 2, which is called as bi-level quantization. In bi-level 
quantization scheme, the sensing threshold can also be considered as the quantized threshold 
to separate the whole value range of test statistics T { into two non-overlapping regions: Region 
and Region I. Correspondingly, there are two quantization levels denoted by q and qi. The 
quantization of T { is formulated by 



where T[ is the quantized test statistics in zth SU. And the report of individual SU is generated 
as following 



In the fusion center, the decision fusion can be treated as that the quantized test statistics is 
combined based on received reports Ri, which is formulated by 




(8) 




(9) 



K 



K 




(10) 



1=1 



1=1 
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where T is the final test statistics in fusion center; k denotes the number of reports R { = 0, 
and k\ denotes the number of reports Ri — 1. Then the final decision is made by 

f H , if T < N (h , if k q + k iqi < N 

Ro(V = \ =\ (ID 

(H u if T > N [H u if k q + k iqi > N 

where the final decision threshold N is belong to the set {/co<?o + k\q± : V/co € {0, 1, • • • , K — 
l}, hi = K — k }. Actually, q and qi can be normalized to f3 = (q — go) / (<?i — qo) = and j3i = 
(<?i — <?o)/(<?i — qo) — 1 separately, regardless of their original values. Thus, T = k x = Y^i Ru 
and we arrive at {k q + kiqi : V/c G {0, 1, • • • , K — 1}, k\ = K — k } := {1, 2, • • • , K}, 

( H , if J2?Ri<N 
Ro=l ■ (12) 

[h u if Ef Ri>N 

When compared (12) with (2), it is evident to see that, from the perspective of quantization 
theory, the conventional decision fusion rule is explained correctly. 

Based on the quantization theory, the cooperative false alarm probability Pf and cooperative 
detection probability P d are expressed by P(T > N\H ) and P(T > N\Hi) respectively. For 
P(T > N\H ), the probability mass function is P(T = N\H ) = P{k q + k iqi = N\H ). Since 
P(koqo + k±qi = N\H ) can be interpreted as the probability of arriving exactly k\ successes in 
K trials, we have 

P(k q + k lQl = N\H ) = kil ( K L ki y P ( T > to\H ) kl (l - P(T > e \H )) K ^ . (13) 

Additionally, because of K = k + hi, k\ = (N — Kq )/(qi — q ). Then, the corresponding 
cumulative distribution function is 



K fK\ 

P(k q + k lQl > N\H y ) J2 ( k ) P ( T ^ e o\H y ) kl (l " P(T > e \H y )) K - h 

k 1 = (N-Kq )/(q 1 ~q ) ^ X ' 



(14) 

where y e {0, 1}. Since q and q 1 can be normalized to /3 — and /3 X = 1, the cumulative 
distribution functions is re-formulated by 

P Y = P(k f3 + kifii > N\H y ) = P(h > N\H y ) 

^ » 1 — * u - \K-ki 



^ ( k ) P ( T >^\H y ) kL (l-P(T>e \H y )y 

k 1= N ^ X ' 



September 15, 2011 



DRAFT 



9 



where Y e {d, /}, y G {0, 1}; Y = d when y — 1, and Y — f when y = 0. When compared 
with (3) (4), (15) proves that the conventional decision fusion rule can be interpreted correctly 
from the perspective of quantization theory. 

B. Performance of CSS with Multi-level Quantization 

From the perspective of quantization theory, the fusion rule with multi-bits report is interpreted 
as the multi-level quantization based hypothesis test. In homogeneity scenario, the test statistics 
Tj is considered to follow a same distribution for different SU. Thus, let T denote the test 
statistics in every SU. Under multi-level quantization, there are I quantization levels with / — 1 
quantization thresholds. Let the I quantization levels and the I — 1 quantization thresholds be 
denoted by {q , q±, - ■ ■ , qi-i} and {e , e±, • • • , q_ 2 } respectively. The / — 1 quantization thresholds 
separate the value range of T into / non-overlapping regions. The probability that T falls in a 
certain region is given as follows, 

P(ei-i <T<e t \H y )= V f(x\H y )dx, for i = 0, 1, • • • , I - 2. (16) 

where f(x\H ) and f(x\H 1 ) denote the pdf's of T under H and Hi respectively; let e_i = — oo 
and e,_i = +oo, hence P(e_i < T < e ) = P(T < e ), P(e,_ 2 <T < e,_i) = P(T > e,_ 2 ). 

The Section III- A shows the case that 1 — 2. Next, we derive the cooperative false alarm 
probability and cooperative detection probability when / = 3. According to quantization theory, 
the two quantized thresholds separate the whole range of T into three regions: Region 0, Region 
1 and Region 2. Correspondingly, there are three quantization levels which are denoted by 
g , ?i, and g 2 respectively. In addition, q , qi, and q 2 can be normalized to f3 — 0, /3i — 
(Qi — 9o)/(?2 — Qo), and /3 2 = 1 respectively. For convenience, f3 , B u and B 2 are directly used 
in following discussion instead of q , qi, and q 2 . When compared with decision fusion rule, the 
value range of N is [0, K\. For the tri-nomial distribution, the related probability mass functions 
are 

P(k B + ki/3i + k 2 /3 2 = N\H y ) = P(ki + k 2 B 2 = N\H y ) 

K \ (IV) 

- — P(T < e \H y ) k0 P(e <T< €i\H y ) k ^P(T > ei\H y ) k * 



k \ki\k 2 \ 

where k , ki, and k 2 denote the number of the value of T falling in Region 0, Region 1, and 
Region 2 respectively, and k + k\ + k 2 = K. 
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Next, we show the conditions that ensure k /3 + kifii + k 2 /3 2 > N holding for any given 
value of N. k /3 + ki/3i + k 2 /3 2 > N is equivalent to hi fa + k 2 > N. We can derive that when 
k 2 > max(0, r^frl) and h > max(0, \^]), 2 h(3i + k 2 > N. Accordingly, the related 
cumulative distribution functions are 

P(k f3 + fcjft + k 2 (3 2 > N\H y ) 

fc2=&2 «i=ai 

where a 2 = max(0, [ ^"l^ 1 ] ), ai = max(0, \ N ^ 2 ])- (18) is a closed form function for coop- 
erative false alarm probability and cooperative detection probability under tri-level quantization. 
Sequentially, we present the cooperative false alarm probability and cooperative detection prob- 
ability under /-level quantization. 

Proposition 1: When / > 2, q < qi < q 2 < ■ ■ ■ < qi-\, f^°° f(x\H )dx < f{x\H\)dx 
for Ve G (— oo, +oo), and ViV G [0, if], the cooperative false alarm probability and cooperative 
detection probability of /-level quantization are given by 

P Y (N) = P(k f3 + kipi + ■■■ + fcz-iA-i > N\H y ) 



(19) 



where 



P(T > e^lHy)*-* . . . P(ti<T<e 2 \H y ) k > 

fei-l=Oi-l fc2=c22 

K-Yj--\ hi 

x £ SSra *<T< e J B .)"P(T<,| J? /, 

fci=ai 



A = (20) 
Qi-i ~ % 



aj = 



max (0, \(N - - ),3=l~l- 

iax (0, [ [N - EtJ +1 hifii -(K- XtJ+i W-i] /(& - j ) , 1 < J < i - 2. 

(21) 

Proof: (19) can be obtained by extending the derivation of functions (15)(18) under / = 2, 3 
to the case / > 3. ■ 
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Mathematically, for N, let the set TZ := {Xl!=o hPi '■ V{ko, fci, • • • , fcj-i}, X^Io ^ = ^} an( ^ 
the maximum number of elements in TZ is Let R denote a matrix including every 

element in TZ in ascending order. From (19)-(21), the following Proposition 2 is derived directly. 

Proposition 2: When iV G [0, K], (19) is the non-increasing function of iV; when N E TZ, 
(19) is the decreasing function of N. 

According to Proposition 2, the cooperative false alarm probability and cooperative detection 
probability can be controlled by adjusting iV in (19), rather than enumerating the permutation of 
{A; , fci,--- , h-i} in (5) for (6) (7). Thus, (19) suits for NP criterion and CDR criterion feasibly. 



C. Study of Quantization based Sensing Performance 

This subsection explains the relation between the combined test statistics and the probability 
mass in detail. And then based on our analysis, we show the computation of sensing performance 
is simplified in Uniform quantization scheme. Let K be the b x I matrix where b = ; 
K(m, n) = k m , n is the element at m-th row and n-th column. Given m, k rn ^, k m ^ • • • , k m ^ means 
a permutation of {k , • • • , where Y^n=o = K. Namely, the rows of K represent every 
realization of user's reports. Since K can be obtained effortlessly, consider K is available. Then 
each realization of combination of quantized test statistics is given by 

Q + = Kq (22) 
where q = [q , q2,-" i Qi~i] T - Actually, Vm, Q + (m) G TZ. Correspondingly, we have 

P+(m) = K\ J] ^ M • (23) 
- LJ - K m,n ! 

n=l v ' 

which is interpreted as the probability mass related to Q + (m). When Q means the matrix that 
the elements of Q + in ascending order, it exists 

P.Q+ = Q (24) 

where is the permutation matrix and it is available. Then, we can obtain 

P.P+ = Py (25) 

Note, for two arbitrary realizations of {fc , fci, • • • , fcj-i}, it is possible that their combined 
quantized test statistics have the same value. In other words, Q may contain elements which 
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have the identical value. For the elements with the same value N in Q, the corresponding index 
is between min(argmin„ |Q(n) — N\) and max(argmin n |Q(n) — N\). Thus, when N <E 1Z, the 
probability mass function for CSS is 

max(arg minn |Q(n)— N\) 

P(k o Po + k 1 1 + --- + k l - 1 l - 1 = N\H y )= Py ^- W 

n=min(arg min„ |Q(n)— N\) 

Consequently, we have the following corollary. 

Corollary 1: When I > 2, g < q 1 < q 2 < ■ ■ ■ < qi-i, f^°° f(x\H )dx < J e +oc f(x\H 1 )dx 
for Ve e (— oo, +oo), and VW e [0, AT], the cooperative false alarm probability and cooperative 
detection probability are given by 

b 

P Y (N) = Py(n) (27) 

n=bo 

where Py(ra) is given by (25), b = { l+I ^ L ), and 

b = min(arg max — — - — -) . (28) 
n I4{n) — iV 

Proof: According to (22)-(26), Corollary 1 is derived evidently. ■ 
Theoretically, (26) is the probability mass function in terms of the combined test statistics 
(k (3 + ki(3i + • • • + ki-ipi-i). From (26) (27), it is clear to see that the value of quantization 
levels affects the quantity and the order of elements in Q, decides the probability mass function 
in turn (26), and determines the final cumulative distribution function (27). Accordingly, when 
the Uniform quantization scheme is employed, the expression and calculation complexity for 
cooperative false alarm probability and cooperative detection probability are simplified and 
reduced respectively. For uniform quantization scheme, the following function is presented for 
performance evaluation. 

Py = Py * Py * • • • * Vy 



s v - 



(29) 

K 

where * denotes the convolution operation; Py is matrix of c-by-1 where c = K(l — 1) + 1. 
Similarly as the derivation of Corollary 1, Corollary 2 is obtained. 

Corollary 2: When I > 2, g < q 1 < q 2 < ■ ■ ■ < qi-i, f(x\H )dx < f(x\Hi)dx 
for Ve G (— oo, +oo), and VW e [0, AT], the cooperative false alarm probability and cooperative 
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detection probability are given by 



c 



P y (iV) = ^P Y (n) 



(30) 



where Py is given by (29) respectively, and 



1 



c = argmax 



R(n) - N ' 



(31) 



n 



Proof: The convolution operation in (29) guarantees that the probability which is the element 
in Py one-on-one correspond to every realization of (k /3 + ki/3i H h fcj_iA-i) which is in 



Additionally, Py can be obtained by employing FFT transform. Let Vq and V\ be the matrices 



of size C x 1 where C = 2 [c/2] . Vi e {1, 2, • • • , K}, V = [pj, T ] T and Pi = [pf , T ] T . 
Py = IDFT(diag(DFT(V y )) K lcxi) = VJ-\diag{WV y ) K l C xi)) 



where DFT(x) and IDFT(x) denote the discrete Fourier transform (DFT) of x and inverse 
discrete fourier transform of x respectively; diag{x) denotes for a diagonal matrix whose diagonal 
entries are same as the elements of vector x in turn; and W is the DFT matrix. Therefore, (32) 
provides an efficient method to obtain Py. 

D. Normal Approximation 

This subsection presents the Normal approximation of cooperative false alarm probability and 
cooperative detection probability. The Normal approximation provides an efficient approach to 
evaluate the performance of multi-level quantization based CSS under given {e , e± : • ■ • , q_ 2 } and 
{g 0) ?i, • • • , Qi-i}- Additionally, the Normal approximation is used for performance optimization 
in next section. Let p^ = [P(T < e \H y ),P(e < T < e^Hy),--- ,P{T > e l _ 2 \H y )} T , (3 = 
[A) j $2, • • • , A-i] T > an d Py = lxxiPy ■ For I > 2, we have mean fj, and standard deviation a 
for state H , and mean /ii and standard deviation cii for state Hi, 



ascending order. Hence, we can arrive at Corollary 2. 



(32) 



W^diagiWVy^lc 



xl 



fly = PyP, 



(33) 




(34) 
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According to Central Limit Theorem (CLT), T denoting the combination of quantized test 
statistics follows that, 

{Af(lM),<ro), H G 
(35) 

where M denotes Normal distribution. 

Considering the continuity correction, we have to refine the approximation by accounting for 
the fact that the multinomial distribution is discrete while the Normal distribution is continuous. 
At first, we take bi-level quantization as an illustration, and then extend to multi-level quantization 
case. As previously explained, the Binomial distribution is treated as the outcomes from bi-level 
quantization. According to the previous analysis, the effective value of N in (15) is positive 
integer, consisting of (3 and Thus, regarding the continuity correction, we have 

P(T > [N/^\/K\H y ) = P(T > \N/l]/K\H y ) = P(T > ([N/l\ + l)/K\H y ) 

« P(T' >([N}+ 0.5)/K\H y ) = P(T >(\N]- 0.5)/K\H y ) (36) 

fx(x\H y )dx 

K K 

where T' G [0, 1]; fj\r(.) denotes the pdf of Normal distribution. Obviously, when K — > +oo, the 
value of (\N]/K — 0.5/ K) can be considered as the continuous positive number for evaluating 
the cooperative false alarm probability and cooperative detection probability. Consequently, we 
arrive at the following Normal approximation for / = 2, 

where 0.5/K is the continuity correction; Q{.) means the Q-function. Sequentially, the following 
corollary is derived, 

Corollary 3: Given I > 2 and q < q 1 < q 2 < ■ ■ ■ < qi-i, for VA^ e [0, if], the Normal 
approximations of cooperative false alarm probability and cooperative detection probability are 
given by 

Py(N) » Q , (38) 



where 



N = arg max (39) 

n R{n)/K — N 
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Proof: At the fusion center, the value of realization of T belongs to set C := {x/K : Vrr e 71} . 
Based on (39), it is observed that P(R(n) — R(n — 1) = (3± — f3 ) has the highest probability 
for Vn > 1. Thus, similarly as the continuity correction 0.5/K in bi-level quantization case, it 
is regarded that the continuity correction denoted by A is as follows for multi-level quantization 
case, 



The derived Proposition 1 provides a feasible bridge between the cooperative sensing per- 
formance and the quantization configurations. On the one hand, (19) can correctly evaluate the 
performance of CSS, when the quantization levels and quantized thresholds are arbitrary values 
under q < q 1 < q 2 < ■■■ < qi-i, on the other hand, the quantization levels and quantized 
thresholds facilitate the understanding of the determination of CSS performance rather than 
likelihood ratio. Hence, we can utilize (19) to reasonably optimize the related quantization 
configuration for high sensing performance. For example, in bi-level quantization based CSS, 
(14) can correctly evaluate the performance of CSS for a certain value of sensing threshold, but 
there exists an optimal sensing threshold which determines the quantization levels for achieving 
the optimal sensing performance [18]. 

A. Minimum Mean Square Error (MMSE) based Quantization Scheme and Optimization 

As a merits of adopting quantization theory to evaluate the CSS, there have been a large 
number of literatures on quantization for direct use. According to [19], the Lloyd's quantization 
scheme can achieve the minimum information loss of quantization. When the given pdf of test 
statistics is g(x), the quantization levels {q , q 1 , • ■ ■ , qi~\} are given by (41) to achieve the MMSE 
caused by quantization. 




2K 



(40) 



Based on (33)-(37)(40), we arrive at Corollary 3. 



IV. Quantization and Optimization 



The mean square error (MSE) under quantization is written by, 




(41) 




i-i 



(42) 
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Substituting (41) into (42) yields 




g(x)dx. 



(43) 



Then, we arrive at the following programming, 




x — 



j e 3 1 xg(x)dx\ 
f^_ i g{x)dx J 



g(x)dx 



(44) 



In CR, there are two pdf's f(x\H ) and f(x\H 1 ) for test statistics T. Hence, g(x) is given by 



where w denotes the cost or weight. The involved cooperative SUs have to quantize the test 
statistics T under H and Hi concurrently so that the fusion center can distinguish the real 
state of primary user with high correct probability. The corresponding visual interpretation is 
that: when treating the outline of the outputs of f(x\H ) and f(x\Hi) in terms of x as the 
waveform of target signal, the objective of quantization is to make the fusion center can recover 
this waveform with minimum distortion. Thus, consider w — 1 in following discussion. (41) (45) 
show a feasible solution to obtain the quantization levels under certain quantization thresholds 

[19] shows the existence and uniqueness of the optimal quantization thresholds {eo, e±, • • • , €1-2} 
from which the optimal quantization levels {go,<?i,"" ,Qi-i} can De obtained directly. Thus, 
we can obtain the optimal sensing configuration including the optimal quantization levels and 
the optimal sensing thresholds. When comparing the exclusive search employed for the optimal 
quantization configuration of Lloyd's scheme in current studies , we show a method to obtain the 
optimum thresholds and quantization levels efficiently. According to [19], when the quantization 
thresholds are optimal, there exists 



g(x) 



/(z|tf )+W(z|#i) 



(45) 



1 + w 



fe + ^ + i)/2,forVje{0,l,-- - ,1-2}. 



(46) 



From (41) (46), we have 




(x - q j - 1 )g(x)dx = 0, for Vj G {2, 3, • • • , I - 2}. 



(47) 



which means, given gj_ 2 and qj-i, the optimal qj makes equation (47) hold. 
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Evidently, J^'^+q^i)^^ ~~ qj-i)d( x )dx is a monotone increasing function of qj. Hence, 
fyiQj) = I(q 3 ~2+q^i)/2( x ~ Qj-i)g( x )dx is quasi-convex function of qj and has the minimal 
value when (pj(qj) > 0. From (41)(46)(47), we have 

xg(x)dx/ / g(x)dx, (48) 

-oo J oo 

gi + (eo) = 2 eo -go + M, (49) 

q+(e ) = argmin 2 (? 2 , qf(e ), ?o"(eo)) 

92 (50) 

st. g+ < g 2 < +oo, 2 (g2,gi f (eo),^(eo)) > 0. 
Let <?iV" ^^l) denote the calculated quantization levels under a given value of e - m 
turn, we can obtain q^_ x given by 

qti( e o) =argmin g£ 2 (e ), <?ts( e o)) 

^ " (51) 

st. g+ 2 < qi-i < +cx3,0,_i(gi_i,g,t 2 (eo),g,t 3 (eo)) > 0. 
Proposition 3: Under > 0, two propositions are given: i) when e is variable, the value 
of the Z-th quantization level monotonically increases as e increasing; ii) when q is variable, 
the value of the Z-th quantization level monotonically increases as q increasing. 

Proof: From (48)-(51), it is obvious that q^_ Y monotonically increases in term of the variable 
eo is increasing, q^_ x monotonically increases in term of the variable q is increasing. ■ 
In addition, when min V; _ 1 = 0, e achieves the optimal value. Summarily, for Ve , obtain 

{Qo( € o),qt(e )r- ■ ,qti( e o)h when 0i-i(g,ti(e o ), g,t 2 (eo), 5it 3 ( e o)) = °> {?o"( e o), qti^o), • • • , 
^"t 1 (e )} are optimal. Thus, the original optimization problem with (21 — 1) variables has been 
transformed to an optimization problem with a variable e . Consequently, by employing the 
bisection method for quasiconvex optimization [21], we can obtain the optimal e at first. 
Sequentiality, by using the optimal e , the other optimal sensing thresholds and the corresponding 
optimal quantization levels can be obtained directly. 

B. Optimization based on Normal Approximation 

Under given quantized thresholds, the following method is proposed to obtain the optimal 
quantization levels by using the results of Normal approximation in Section III-D. From (38), 
we have 

P d = Q (fo) (52) 
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where 

j = Q~ 1 (P/)^o + /io-/ii (53) 

The smaller value of f a gains the higher detection probability. Letting q = [go,<?2,-" ,1i-i] T 
replace (3 in (33)-(34) yields 

h, = pJq> (54) 



q^(I - P y )T(I - Py )q (55) 



= i, K 



Substituting (54)(55) into (53) yields 



f ( g- 1 (i J /)go-(/ii-/io) Q-'Wff- Pq) t (i - Po)q/^ - q T (pi - go) 

*i v /q^(I-P i r(I-P 1 )q/K 

(56) 

Thus, the following programming is obtained, 

min / (q) st. q ^ (57) 
q 

which can be efficiently solved by using the method introduced in [16]. Consequently, given 
sensing thresholds, the optimal quantization levels can be obtained from our programming 
efficiently. 

V. Examples 

For clarity of expression, let LO denote the results of CSS in Lloyd's quantization scheme with 
the proposed optimization method; NA means the results of CSS from Normal approximation; AF 
represents the sensing results when T is analog forwarded to fusion center without quantization; 
LRT expresses the results under LRT test and the quantized thresholds are same as the LO's; 
and let NO denote the results of Normal- Approximation-based-optimization at the case that the 
quantized thresholds are same as the LO's. 

A. Normal Data 

Let T follow the Normal distribution under both states H and Hi, 

1 (x- fi T ,o) 2 , co . 
f(x\H ) = — exp( ^— , (58) 



\'2rraTo °' 



2 

Tfl 
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/(Z|#l) 



1 



exp( 



(x - /iT,l) 2 



(59) 





where /j, Tfi , pL T ,i, ct.o* and cr T ,i are the mean and standard deviation under H and Hi respectively. 
Fig. 2-Fig. 7 show the results of cooperative sensing under quantization with fi T:0 = 0, fi T ,i = 1, 
<tt,o = 1-5, and cr Tjl = 2. 3 It is apparent to obtain the following conclusions: i) the LO case 
can achieve the high sensing performance when compared with the AF case; ii), according to 
Fig. 2, Fig. 4, and Fig. 6, the Normal approximation provides good result for analyzing the 
CSS's performance under quantization. And the results of Normal approximation become more 
accurate as the quantity of quantization levels or the number of cooperative users increasing; Hi) 
the NO case can achieve comparable sensing results to LRT case, according to Fig. 3, Fig. 5, 
and Fig. 7. 

B. Rayleigh Data 

Consider T follows Normal distribution under H , and Rayleigh distribution under Hi. The 
related pdf under Hi is given by 



Fig. 8-Fig. 11 indicate the results of CSS under quantization with fi Tfi = 1, cr Tj0 = 0.5, and 
<7t,i = 1. It is evident that we can obtain the conclusions same as the Normal data case's, even 
the Rayleigh data is employed. Hence, it proves that our derived methods and functions can suit 
for different distributions of test statistics. 



In this paper, the performance of CSS is studied considering the impact of quantization. 
We derive the closed form functions to calculate the cooperative false alarm probability and 
cooperative detection probability, which ensure the performance of CSS with limited quantization 
levels is no longer intractable, especially, for optimization. Additionally, the related optimization 
methods are provided to improve the performance of CSS. 

3 In practice, the values of /lit.o and fir,i can be normalized to and 1 respectively. Here, we prefer <jt,o = 1.5 and ctt,i = 2 
randomly with considering that <tt,o < ot,\ in CR scenario. For example, the CSS employs energy detector. 




(60) 



VI. Conclusion 
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Fig. 2. The ROC line with K — 3, I = 3. Marks denote the CSS's performance under the different discrete combined test 
statistics. 
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Fig. 3. The sensing performance under different N, with K = 3, I = 3. Marks denote the CSS's performance under the 
different discrete combined test statistics. 




Fig. 4. The ROC line with K = 5, I = 3. Marks denote the CSS's performance under the different discrete combined test 
statistics. 
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Fig. 5. The sensing performance under different N, with K = 5, I = 3. Marks denote the CSS's performance under the 
different discrete combined test statistics. 




Fig. 6. The ROC line with K — 3, I = 4. Marks denote the CSS's performance under the different discrete combined test 
statistics. 
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Fig. 7. The sensing performance under different N, with K = 3, I = 4. Marks denote the CSS's performance under the 
different discrete combined test statistics. 




Fig. 8. The ROC line with K = 3, I = 3. Marks denote the CSS's performance under the different discrete combined test 
statistics. 
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Fig. 9. The sensing performance under different N, with K = 3, I = 3. Marks denote the CSS's performance under the 
different discrete combined test statistics. 




Fig. 10. The ROC line with K = 3, I = 4. Marks denote the CSS's performance under the different discrete combined test 
statistics. 
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Fig. 11. The sensing performance under different N, with K = 3, I = 4. Marks denote the CSS's performance under the 
different discrete combined test statistics. 
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